arXivil504.06693vl [quant-ph] 25 Apr 2015 


Controlled Electromagnetically Induced Transparency and Fano Resonances in Hybrid 

BEC-Optomechanics 

Kashif Ammar Yasir* and Wu-Ming Liu^ 

Beijing National Laboratory for Condensed Matter Physies, 

Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China. 

(Dated: April 28, 2015) 

We investigate the controllability of electromagnetically induced transparency (EIT) and Fano 
resonances in hybrid optomechanical system which is composed of cigar-shaped Bose-Einstein con¬ 
densate (BEC) trapped inside high-finesse Fabry-Perot cavity driven by a single mode optical held 
along the cavity axis and a transverse pump held. Here, transverse optical held is used to control 
the phenomenon of EIT in the output probe laser held. The output probe laser held can efficiently 
be amplihed or attenuated depending on the strength of transverse optical held. Furthermore, we 
demonstrate the existence of Fano resonances in the output held spectra and discuss the controlled 
behavior of Fano resonances using transverse optical held. To observe this phenomena in laboratory, 
we suggest a certain set of experimental parameters. 
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I. INTRODUCTION 

Over a few years, Cavity-optomechanics, a rapidly 
developing area of research, has made a remarkable 
progress. A stunning manifestation of optomechanical 
phenomena is in exploiting the mechanical effects of light 
to couple the optical degree of freedom with mechanical 
degree of freedom [1]. In this regard, a milestone was 
achieved by the demonstration of optomechanics when 
other physical objects, most notably cold atoms or Bose- 
Einstein condensate (BEC), were trapped inside cavity- 
optomechanics [2, 3]. Through several investigations, dif¬ 
ferent aspects of BEC have made spectacular contribu¬ 
tions in understanding the complex systems [5-8]. In op¬ 
tomechanical systems coupling is obtained by radiation 
pressure [9, 10] and indirectly via quantum dot [11] and 
ions [12]. Optomechanics helps mechanical effects of light 
to cool the movable mirror to its quantum mechanical 
ground state [13-17] and provides a platform to study 
strong coupling effects in hybrid systems [18-20]. Re¬ 
cent magnificent discussions and simulations on bistable 
behavior of BEC-optomechanical system [21], high fi¬ 
delity state transfer [22, 23], entanglement in cavity- 
optomechanics [24-28], dynamical localization in field 
of cavity-optomechanics [29] and the coupled arrays of 
micro-cavities [30] provide clear understanding for cavity- 
optomechanics. 

Recently, electromagnetically induced transparency 
(EIT), a phenomenon of direct manifestation of quantum 
coherence [31, 32], has been extensively investigated and 
has provided a lot of remarkable applications [33, 34]. 
In atomic system, the EIT occurs due to quantum in¬ 
terference effects induced by coherently driving atomic 
wavepacket with an external pump laser field [35]. The 
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EIT effect has been theoretically investigated in optome¬ 
chanical system [36, 37] and later experimentally verified 
in both optical [38, 39] and microwave [40] domains. The 
Fano resonances [41] have played an important role in 
understanding the photo-electron spectra [42] in atomic 
physics and have made a magnificent contribution in the 
latest field of plasmonics [43]. Fano resonances have also 
been investigated in hybrid cavity-optomechanics by us¬ 
ing different configuration [44, 45]. 

In this paper, we investigate the controlled behavior 
of electromagnetically induced transparency (EIT) and 
Fano Resonances in hybrid BEC-optomechanical system 
which is composed of cigar-shaped Bose-Einstein conden¬ 
sate (BEC) trapped inside high-finesse Fabry-Perot cav¬ 
ity, with one fixed mirror and one moving-end mirror 
with maximum amplitude driven by a single mode 
optical field r] along the cavity axis and a transverse 
pump field r ]±. Transverse optical field is used to control 
the phenomenon of electromagnetically induced trans¬ 
parency (EIT) in the output probe laser field. By observ¬ 
ing output probe field spectra, we show that the probe 
laser field can efficiently be amplified or attenuated de¬ 
pending on the strength of transverse optical field r]_\_. 
Furthermore, we demonstrate the existence of Fano res¬ 
onances by solving output field spectra and discuss the 
controlled behavior of Fano resonances using transverse 
optical field. To observe this phenomena in laboratory, 
we have suggested a certain set of experimental parame¬ 
ters. 

In this paper, the mathematical model of the hybrid 
BEC-optomechanical system is presented in section II. In 
section III, we derive the Langevin equation to introduce 
the dissipation effects in the intra-cavity field, damping 
associated with moving-end mirror and depletion of BEC 
in the system via standard quantum noise operators. In 
section IV, we discuss EIT in output probe field of the 
hybrid cavity-optomechanics and explain controlled be¬ 
havior of EIT with transverse optical field. In section V, 
we demonstrate the existence of fano resonances and de- 
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FIG. 1: (Color online) Cigar-shaped Bose-Einstein conden¬ 
sate (BEC) is trapped inside a Eabry-Perot cavity of length L 
with a fixed mirror and a moving-end mirror, with frequency 
Um and having maximum amplitude of go, driven by a single 
mode optical field with frequency cUp. Intracavity field, with 
frequency c^c, interacts with BEC and due to photon recoil, 
the momentum side modes are generated in the matter wave. 
To observe electromagnetically induced transparency (EIT) 
in output field spectra, we use external probe field with fre¬ 
quency LUpp. Ep corresponds to the power of external probe 
field. We use another transverse field, with intensity 77 ^, to 
control the output spectral dynamics of optomechanical sys¬ 
tem. 


scribe their controllability by using transverse optical 
field. The results are summarized in section- VI. 


II. BEC-OPTOMECHANICS 


We consider a Eabry-Perot cavity of length L with a 
fixed mirror and a moving-end mirror driven by a single 
mode optical field with frequency as shown in Fig. 1. 
Cigar-Shaped BEC, with N-two level atoms, is trapped 
inside optical cavity [2, 46, 47]. Moreover, optomechani¬ 
cal system consists of a transverse optical field, with in¬ 
tensity r]± and frequency 00 which interacts perpen¬ 
dicularly with BEC [48-52]. Counter propagating field 
inside the cavity forms a one-dimensional optical lattice. 
The moving-end mirror has harmonic oscillations with 
frequency oOm having a maximum amplitude go and ex¬ 
hibits Brownian motion in the absence of coupling with 
radiation pressure. 

The complete Hamiltonian of the system consists of 
three parts, 

H = + ( 1 ) 

where, Hrn is related to the intra-cavity field and its cou¬ 
pling to the moving-end mirror, Ha describes the BEC 
and its coupling with intra-cavity field while, Ht ac¬ 
counts for noises and damping associated with the sys¬ 
tem. The Hamiltonian Hrn is given as [53], 

Hra = + q^) - gmhc^cq-ihri{c-o') 

- ihEp{ce^^^ ( 2 ) 


where, first term corresponds to the energy of the intra¬ 
cavity field, Ac = ujc — oOp is detuning, ujc is intra-cavity 
field frequency and (c) are creation (annihilation) op¬ 
erators for intra-cavity field interacting with mirror and 
condensates with commutation relation [c, c^] = 1 . Sec¬ 
ond term accounts for the energy of moving-end mir¬ 
ror. Here g and p are representing dimensionless position 
and momentum operators respectively, for moving-end 
mirror with commutation relation [g,p] = i, which re¬ 
veals the value of the scaled Planck’s constant, h = 1. 
As intra-cavity field is coupled with moving-end mir¬ 
ror through radiation pressure so third term represents 
coupled energy of moving-end mirror with field. Here 
Pm = V^{uJclL)xo is the coupling strength and xq = 
^yh/2mujrn, is zero point motion of mechanical mirror 
having mass m. Second last term gives relation be¬ 
tween the intra-cavity field and the external pump field 
1 77 1 = PinK/hjjp^ where, P is the input laser power 
and K is cavity decay rate associated with out-going 
modes. Last term is associated with external probe field 
and Ep is related to the power of the probe field as 
\Ep\ = Pphi/hjujp. Ap = ujp — uji is the detuning of 
external probe laser field with external pump field. 

The Hamiltonian for BEC and its coupling with intra¬ 
cavity field (i7a), in strong detuning regime and in the 
rotating frame at external field frequency is derived by 
considering quantized motion of atoms along with the 
cavity axis. We assume that BEC is dilute enough and 
many body interaction effects are ignored[49, 51]. Here, 

Ha = [ i^^{x)[- P HUoc^ccos^ kx 

J 2madx^ 

+ hri^cos{kx){c^c)\ip{x)dx^ (3) 

is bosonic annihilation (creation) operator and 
^0 = go is the vacuum Rabi frequency, is 

far-off detuning between field frequency and atomic tran¬ 
sition frequency. [Note: in equation 3, we have not con¬ 
sidered for now the effects of harmonic trap causing the 
confinement of atoms inside the cavity.] Furthermore, 
rria is mass of an atom, ujq is atomic transition frequency 
and k = ujp/c is the wave number. 77 ^ = goQp/Aa is 
the coupling of BEC with transverse field and represents 
maximum scattering and ftp is the Rabi frequency of the 
transverse pump field. Due to field interaction with BEC, 
photon recoil takes place that generates symmetric mo¬ 
mentum E2lhk side modes, where, I is an integer. In the 
weak field approximation, we consider low photon cou¬ 
pling. Therefore, only lowest order perturbation of the 
wave function will survive and higher order perturbation 
will be ignored. So, 7 /; is defined depending upon these 
Qth^ fst 2'^^ modes [21] as, 


ip{x) 


1 

7z 


[boV2cos{kx)bi P V2cos{2kx)b2], (4) 


here, bo^bi and 62 are annihilation operators for 0 ^^ , 1 ^^ 
and 2'^^ modes respectively. By using ip{x) defined above 
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in Hamiltonian Ha^ we write the Hamiltonian governing 
the field-condensate interaction as, 

Ha = uJrbibi -f 4:UJrb2^b2 + c[V2{bo^b2 + b2^bo) 

+ 2N + +c}[V2{bohi 

+ bi bo) + {bi 62 + 62 61 )]. (5) 

The sum of particles in all momentum side modes is, 

+ 61^61 + 62^62 = N, where, N is the total num¬ 
ber of bosonic particles. As population in 0 ^^ mode is 
much larger than the population in and 2'^^ order 
side modes, therefore, we can comparatively ignore the 
population in and 2 ^^ order side modes and can write 

60 60 — or 60 and 60 ^ \/N. This is possible when 
side modes are weak enough to be ignored. Moreover, 
for Bogoliubov mode expansion, we consider small in¬ 
teraction of intra-cavity field with BEC so, that atomic 
mode can perform motion analog to the moving-end mir¬ 
ror of the system. Under these assumptions, we recover 
the cavity-optomechanics-like Hamiltonian discussed in 
Ref. [2] and [21] and given as, 

Ha = c^c+^{P^+Q^)+gahc^cQ+fi7]effQ{c^+c). 

(6) 

First term accommodates the potential energy for con¬ 
densate in intra-cavity field. We assume large atom-field 
detuning A^, so that, the excited atomic levels can be 
adiabatically eliminated. Second term describes the mo¬ 
tion atomic momentum side modes exited by radiation 
pressure. It can be observed that the atomic side modes 
are analogous of a mirror whose motion is driven by radi¬ 
ation pressure. P = ^(62 — ^) and Q = ^(62 + 62 ) are 
dimensionless momentum and position operators for such 
atomic mirror with canonical relation, [Q, P] = i and Cl = 
4uJr = 2hk^/ma, is recoil frequency of an atom. Third 
term in equ.( 6 ) describes coupled energy of field and con¬ 
densate with coupling strength Qa = 
where, rrihec = {L‘^NUoOJr) is the side mode mass 
of condensate. The last term accounts for the coupling 
of BEC with transverse field and rjef f = is trans¬ 

verse coupling strength. From equ. 6 , it can be noted that 
in the absence of transverse optical field, when there is 
no excitation for cos{kx)^ we recover the same expression 
for atomic mode of the cavity-optomechanics as in Ref. [2] 
and [ 21 ]. 


III. LANGEVIN EQUATIONS 

The Hamiltonian Ht describes the effects of dissipa¬ 
tion in the intra-cavity field, damping of moving-end mir¬ 
ror and depletion of BEC in the system via standard 
quantum noise operators [54]. The total Hamiltonian 


H leads to develop coupled quantum Langevin equations 
for optical, mechanical (moving-end mirror) and atomic 
(BEC) degrees of freedom. 


dc 

dt 

dp 

dt 

dq 

dt 

dt 

dQ 

dt 


= c = (iAo + igmq - WaQ - i^)c + irjeffQ 


Tj + EpC 

( 7 ) 

p = -Umq + gmC^C - JrnP + / b , 

(8) 

q = UJmP, 

( 9 ) 

P = - 4 uJrQ - gaC^C - -frP + /im, 

(10) 

Q — 4 uJj,P 7rQ + /2m- 

(11) 


A = Ac — NUo/2 is the effective detuning of the system 
and din is Markovian input noise associated with intra¬ 
cavity field. The term 7^1 describes mechanical energy 
decay rate of the moving-end mirror and fs is Brownian 
noise operator associated with the motion of moving-end 
mirror [55] . The term 7 ^ represents damping of BEC due 
to harmonic trapping potential which affects momentum 
side modes while, firn and / 2 m are the associated noise 
operators assumed to be Markovian. We consider posi¬ 
tions and momenta as classical variables. To write the 
steady state values of of the operator, we assume optical 
field decay at its fastest rate so that the time derivative 
can be set to zero in equation (7). The static solutions 
are given as. 


Cs 


Qs 

Ps 

Qs 

Ps 


TJ “h iTJgjj^Q 

K + i{Aa - gmQ + 9aQ)' 

gmC^C 

UJfn 

0, 


-gqC^C 

4uJr[^ — ^r/4uJr] ’ 



( 12 ) 

(13) 

(14) 

(15) 

(16) 


where c^, Qs and Qs represent the steady-state solution 
of intra-cavity field, the mechanical mirror displacement, 
and the position of the BEC mode, respectively. To ob¬ 
serve output field spectra, we deal with the mean re¬ 
sponse of the system to the probe field in the presence 
of the coupling field and transverse field. First, we lin¬ 
earized quantum Langevin equations by inserting ansatz 
c{t) = Cs + 6c{t), q{t) = gs+ Sq{t), p{t) = Ps + 5p{t), 
Q\t) = Qs + SQ{t) and P\t) = Pg + SP{t) in Langevin 
equations and taking care of only first-order terms in 
fluctuating operators 6 c(t), Sq{t), bp(t)^ SQ{t) and SP{t). 
The linearized quantum Langevin equations are obtained 
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as, 


6c{t) = -{hv ^ iA)6c{t) ^ GmSq{t) - GaSQ{t) 

'^Veff^Q{t) + EpC (17) 

^ct(^) = - Ga^^Q(^) 

- irjeffSQ{t) + + V^Cin, (18) 

6q{t) = uJmSp{t), (19) 

Sp(t) = -UJm^q{t) + Gm{^c{t) + (^c’*‘(t)) 

- 7mSp{t)EfB, (20) 

SQ{t) = UJrSP{t) + f2m, (21) 

Sp{t) = -UJrn^Q{t) E Ga{Sc{t) E 6p{t)) 

— Jr^P{t) + /im: (22) 

where, A = Aq — gmq E 9aQ is the effective detuning 
of the system and Gm = 9m\cs\, Ga = ga\cs\ are the 
effective coupling of optical field with the moving-end 
mirror and the condensate mode, respectively. To solve 
mean value equation of the system, we write expectation 
value of operators in form O = here O 

is a generic operator. < 6c{t) >, < Sq{t) >, < Sp{t) >, 
< SQ{t) > and < SP{t) > are the expectation values 
corresponding to fiuctuating operators Sc{t), ^q{t)^ ^p(t), 
6Q{t) and SP{t), respectively. 

To solve linearized quantum Langevin equations, we 
assume that the coupling of external pump field g is much 
stronger than the coupling of external probe field Ep. Un¬ 
der this assumption, the solutions of linearized Langevin 
equations can be approximated to the first order external 
probe field by ignoring all higher order terms of Ep. So, 
the solution for intra-cavity probe field is now given as. 


En 


c+ = 


- [Y(Ap) + (k -h i(A Y Ap) 


X{Ap)Y{Ap)^ 

Y X{Ap)){K-i{A-Ap)YX{Ap))l 
Kj — i(A + Ap) -T A*(Ap) 


C- = Er 


y*(Ap) 


where. 


X(Ap) = - 


_ ^a^'f' '^^aVef f^r 


UJ^ - AIy i^rAp 


G‘^idm 


- AlYijmAp' 

Y{Ap) = -A‘^ - iP-2iBApY Al 

2ujj-G -|- iAp) 2ujrpAgQjjGa 
ijr^p - AjYuJ^ 

2KjUJmG 


ijmAp - A2 


(23) 

(24) 


(25) 


(26) 


Eq.23 and Eq.24 clearly describe the dependence of 
output field spectra on the coupling of different degrees of 
freedom. In particular, we can observe the rule of trans¬ 
verse optical field coupling with BEG mode in output 
field. Eurthermore, to investigate the EIT-like behavior. 


we write output field spectra by using input-output rela¬ 
tion Cout = V^c — Cin [54] , where and Cout represent 
input and output field, respectively. Moreover, we ignore 
quantum noises associated with Cout and Cin as discussed 
earlier. The out-going optical field can be expressed as, 

< Cout > = Co + c+Epe-^^-^ + (27) 

By using above relation and input-output field relation, 
we describe the components of output field spectra at 
probe field frequency and Stokes frequency as. 


c+ = -h 1, 

c_ = v^c_, (28) 


respectively. In the absence of optical coupling with 
moving-end mirror and BEG mode, the output field spec¬ 
tra at probe frequency and stokes frequency is given as. 


2k 

K “b ^(A — Ap) 
c+ = 0. 


(29) 

(30) 


IV. CONTROLLABLE EIT IN OUTPUT FIELD 

The hybrid BEG-optomechanical system shown in 
Eig.l is simultaneously driven by external pump field 
with frequency uOp and probe field with frequency cjpp, 
generating radiation pressure force which oscillates at 
frequency difference Ap = ujpp — ujp. When this resul¬ 
tant force oscillates with frequency close to the frequency 
of mechanical mode (moving-end mirror) uum or atomic 
mode (BEG) it gives rise to the Stokes and anti- 
Stokes scatterings of light from the strong intra-cavity 
standing field. But the Stokes scattering is strongly 
suppressed because conventionally, optomechanical sys¬ 
tems are operated in resolved-sideband regime k « ujmi 
which is off-resonant with Stokes scattering and so only 
anti-Stokes scattering survives inside the cavity. There¬ 
fore, due to the presence of probe field and pump field 
inside the system, electromagnetically induced trans¬ 
parency (EIT) like behavior appears in the output field 
spectra. 

To make this study of tunable EIT and Eano res¬ 
onances in hybrid BEG-Optomechanics experimentally 
feasible, we choose a regime of particular parameters very 
close to the recent experimental studies [1, 2, 47]. In ad¬ 
dition, we consider the parameters such that the system 
remain in stable regime as discuss in Ref. [50, 51]. We 
consider N = 2.3 x 10^ atoms trapped inside Eabry- 
Perot cavity with length L = 1.25 x 10“^, driven by sin¬ 
gle mode external field with power Pin = 0.0164mIU, 
frequency ujp = 3.8 x 27r x 10^^Hz and wavelength 
Xp = 780nm. The intra-cavity optical mode oscillates 
with frequency ooc = 15.3 x 27r x with intra-cavity 

decay rate k = 1.3 x 2'KkHz. The vacuum Rabi frequency 
of the system is considered as, Uq = 3.1 x 27rMHz. 
Intra-cavity field produces recoil of = 3.8 x 27rkHz 









5 




FIG. 2 : (Color online) The absorption and dispersion quadra¬ 
tures of the output field spectra C- are demonstrated as 
a function of probe detuning /\plujm and in the absence 
of of transverse optical field coupling rjeff/t^ — 0. (a) de¬ 

scribe the behavior of single-EIT in real quadrature of out¬ 
put probe field with various coupling strengths Gal^m — 0 , 
Gml^m = 0 (blue curve), Gal^m = 0.02, Gml^m = 0 (red 
curve), and Gal^m — 0.03, Gml^m — 0 (green curve), (b) 
contains similar behavior in imaginary quadrature of output 
probe field with different couplings Gal^Um = 0 , Gml^m — 0 
(blue curve), Gal^m — 0.02, Gml^m — 0 (red curve), and 
Gaj^m — 0.03, Graj^m — 0 (grccu curvc). Thc effective de¬ 
tuning of the system is A = 0.52 x 2'kMHz and the frequency 
of moving-end mirror is considered as, ujm — 1-02 x 2'kMHz. 
The atomic mode and mechanical mode damping are 7 r = 
0.21 X 2KkHz and 7 m = 1.1 x 2KkHz, respectively. The cav¬ 
ity optical decay is /c = 1.3 x 2KkHz. 


in atomic mode trapped inside cavity with damping rate 
7 ^ = 0.21 X 2KkHz. The moving-end mirror of cavity 
should be a perfect reflector and oscillates with frequency 
uJm = 1.02 X 2kMHz with damping 7 ^, = 1.1 x 2KkHz. 
From given parameters, one can observe that the system 
is in resolved-sideband regime because ujrn >> /^, this 
condition is also referred to good-cavity limit. 

The real {Re[C-\) and imaginary (/m[C_]) quadra¬ 
tures of out-going probe field as a function of probe- 
cavity detuning are discussed in Fig2, in the absence 
of transverse field coupling rjeff/n = 0. Here, Re[C-\ 
and Im[C-\ accounts for in-phase and out-phase, respec¬ 
tively, quadratures of output field spectra and also re¬ 
ferred to the absorption and dispersion behavior of out¬ 
going optical mode. Fig. 2 (a) and Fig. 2(b) demonstrate 


the single-EIT behavior of such absorption and disper¬ 
sion quadratures, respectively, of output field spectra for 
different coupling strengths Gaj^m = 0 , Gml^m = 0 
(blue curve), Gal^m = 0 . 02 , Gml^m = 0 (red curve), 
and Gaj^m — 0.03, Gml^m = 0 (green curve). To ob¬ 
serve single-EIT behavior, we have considered the case 
when system is only coupled to the intra-cavity atomic 
mode or BEG and the coupling of moving-end mirror 
with intra-cavity optical mode is zero. The EIT behavior 
in cavity-optmechanics with coupled mirror have already 
been discussed in previous works like [36, 38, 39, 44]. We 
consider that the optomechanical system is being oper¬ 
ated in strong coupling regime which means intra-cavity 
optical mode is strongly coupled to atomic mode of the 
system. When collective density excitation of atomic 
mode of cavity became resonant with intra-cavity opti¬ 
cal mode, it cause strong coupling between atomic mode 
and system. It is only possible when the strength of 
coupling between single atom and single photon of the 
cavity go = 10.9 x 2kMHz is larger than both decay 
rate of the atomic excited state 7 ^ = 0.21 x 2KkHz 
as well as intra-cavity field decay k. = 1.3 x 2KkHz 
(^0 >> 7 a 7 ^)[ 2 , 39]. We can also note from mathe¬ 
matical expression of atomic coupling that the strength 
of atomic mode coupling is directly proportional to the 
vacuum Rabi frequency {Ga ocUo = go/Aa). 

In Fig. 2, one can easily observe there are no signatures 
of EIT in absorption and dispersion spectra of output 
field (blue curves) when intra-cavity field is not coupled 
with mechanical mode and condensate mode {Gal^m — 
0, Gml^m = 0). While in red curve, the single-EIT 
window appears in output probe field due to intra¬ 
cavity optical mode coupling with atomic mode (BEG) 
{GaloJm = 0 . 02 ), however, coupling of intra-cavity field 
with mechanical mode is kept zero Gml^m = 0. Green 
curve shows the enhancement in single-EIT window in 
output probe field by increasing the coupling strength to 
Gal^m = 0.03. These results clearly prove the existence 
of single-EIT window in output probe field when intra¬ 
cavity optical mode is only coupled to atomic mode of 
cavity-optomechanics. 

Fig. 3 describes the EIT behavior in the output field 
spectra of the optomechanical system in presence of 
probe detuning Ap/ujm and transverse optical field 
77e////^. Fig. 3(a) and Fig. 3(d) represent absorption (real) 
and dispersion (imaginary) quadratures, respectively, of 
output probe field in the absence of transverse field 
rjeff/K = 0 with various coupling strengths Gaj^ra — O 7 
GjYifuJm — 0 (blue curve), Gaf^m — 0.08, GfYifujjyi — 0.05 
(red curve), and Gal^m — 0 . 1 , Gml^m = 0.08 (green 
curve). Results show that there are no signs of EIT- 
like behavior in absorption and dispersion spectra of 
output field (blue curves) when system is isolated form 
mechanical mode and condensate mode [Gat^m — 0 , 
Gml^m = 0). On the other hand, in red curve, two 
EIT windows appear in output probe field because the 
optical mode of the system is now coupled to both me¬ 
chanical mode (moving-end mirror) {Gmt^m = 0.05) as 
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FIG. 3: (Color online) The real and imaginary quadratures of the output probe field C- are presented as a function of probe 
detuning Ap/curn and transverse optical field r/effli^- (a) and (d) accommodate double-FIT behavior in real and imaginary 
quadratures, respectively, in the absence of transverse field ?7e///^ = 0 and with various coupling strengths Gajojrn — 0, 
Gml^m = 0 (blue curve), Gaj^m — 0.08, Gml^m = 0.05 (red curve), and Gal^m — 0.1, Gml^m = 0.08 (green curve), (b) 
and (e) demonstrate single-EIT windows with coupling strengths Gajojm — 0.03, Gml^Jm = 0, as a function of transverse field 
strengths r\efflK — 0 (blue curve), r]efflK, = 0.02 (red curve) and r\efflK — 0.03 (green curve). Similarly, (c) and (f) show 
double-EIT behavior of output field spectra with coupling Gal^Jm — 0.1, Gml^Jm — 0.08, as a function of transverse field 
strengths r]efflK = 0 (blue curve), r]eff Ik, = 0.02 (red curve) and r]eff Ik = 0.03 (green curve). The remaining parameters are 
same as in Fig. 2. 


well as to the condensate mode with coupling strength 
Gal^Jm = 0.08. Such behavior is also known as double- 
EIT response of output field [44, 45]. When system is 
coupled to atomic mode and mechanical mode at the 
same time and optical mode of the system become reso¬ 
nant to both these modes, it gives rise to anti-Stokes scat¬ 
tering inside the system causing appearance of another 
EIT window in output field. Green curves in Fig. 3 (a) 
and (d) demonstrate similar behavior double-EIT when 
the coupling strengths are increased to Gal^m — 
Gml^m = 0.08. We observe that the quadratures of 
double-EIT behavior are increased by increasing coupling 
strengths. Given results in Fig. 3(a) and (d) show such 
double-EIT behavior in output field when optomechani¬ 
cal system is coupled to moving-end mirror of the system 
and BEG trapped inside the system. 

Fig. 3(b) and Fig3(e) demonstrate single-EIT behav¬ 
ior in absorption and dispersion quadratures, respec¬ 
tively, of output probe field under the influence of trans¬ 
verse optical field rjef /Ik when intra-cavity optical mode 
is coupled to condensate mode with coupling strength 
Gal^m = 0.03 while the coupling of optical mode with 
moving-end mirror is zero Gml^m = 0. [Note: we cannot 
observe transverse optical field effects on single-EIT when 
intra-cavity optical degree of freedom is only coupled to 
the moving-end mirror, as shown in single-EIT results 


in previous works like [36, 38, 39, 44], because trans¬ 
verse optical field is only interacting with BEG trapped 
inside the cavity. Therefore, we only consider conden¬ 
sate mode coupling while studying effects of transverse 
field on single-EIT.] Blue curves show single-EIT win¬ 
dows in output probe field in the absence of transverse 
optical field rjeffli^ = 0. On the other hand, red and 
green curves demonstrate the effects of transverse field 
strengths rjefflf^ = 0-02 and rjeffiK = 0.03, respectively, 
on the single-EIT behavior. When transverse field pho¬ 
ton interacts with atomic mode of the system, it gives rise 
to the total photon number n inside the cavity by scat¬ 
tering transverse photons into the system, which leads to 
another nonlinear contribution to the anti-Stokes scatter¬ 
ings and enhance the EIT behavior in output field. We 
can observe such effects of transverse field in the results 
that the strength of single-EIT is efficiently amplified by 
increasing the strength of transverse optical field. 

Similarly, Fig. 3(c) and Fig3(f) represent double-EIT 
behavior in absorption and dispersion quadratures re¬ 
spectively, of output probe field as a function of trans¬ 
verse optical field r]ef f I^ when intra-cavity optical mode 
is coupled to both condensate mode with coupling 
strength Gal^m — 0-1 and to the moving-end mirror 
is Gml^m = 0.08. Blue curves in both these figures de¬ 
scribe double-EIT behavior in the absence of transverse 
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FIG. 4: (Color online) Single-Fano resonances are shown in 
the absorption (real) and dispersion (imaginary) profile of 
output probe field in the absence of transverse field coupling, 
as a function of normalized probe field detuning Aplum and 
normalized effective detuning of the system Ajum with cou¬ 
pling Gal^m — 0.03 and Gal^m — 0. (a) demonstrates ab¬ 
sorption spectra of output field with different values of ef¬ 
fective detuning Ajum = 0.84 (blue curve), Ajum = 0.87 
(red curve) and Ajujm = 0.90 (green curve). Similarly, (b) 
shows dispersion behavior of output field with effective detun¬ 
ing Ajum = 0.99 (blue curve), Ajum = 1.02 (red curve) and 
Ajujyn — 1-06 (green curve). All other parameters are same 
as in Fig. 2. 



FIG. 5: (Color online) Double-Fano resonances are presented 
in the absorption (real) and dispersion (imaginary) profile 
of output probe field in the absence of transverse field cou¬ 
pling, as a function of normalized probe field detuning Apjujra 
and normalized effective detuning of the system Ajum with 
coupling Gal^m — 0.1 and Gaj^m = 0.08. (a) shows real 
quadrature of output field under different effective detuning 
strengths Ajum = 0.83 (blue curve), Ajujm — 0.85 (red 
curve) and Ajum — 0.97 (green curve). On the other hand, 
(b) presents imaginary quadrature of output field with effec¬ 
tive detuning values A/cJrn = 1-05 (blue curve), A/c^m = 1.1 
(red curve) and A/ujm = 1.15 (green curve). Remaining pa¬ 
rameters are same as in Fig. 2. 


field Tjeff/ti = 0. Besides, red and green curves represent 
double-EIT with transverse field strengths = 0.02 

and = 0.03, respectively. We can observe, like 

single-EIT results, double-EIT windows are enhanced by 
increasing the transverse optical coupling. Therefore, in 
accordance to these results, we can confidently state that 
by increasing transverse optical field coupling, we can 
enhance the phenomenon of electromagnetically induced 
transparency (EIT) in hybrid BEC-optomechanics. 


V. TUNABLE FANO RESONANCES 

The formation of Eano resonance in the output opti¬ 
cal mode of hybrid optomechanical system is a fascinat¬ 
ing phenomenon caused by quantum mechanical interac¬ 
tion between different degrees of freedom of the system 
[44, 45] . The constructive and destructive quantum inter¬ 
ferences among narrow discrete intra-cavity optical reso¬ 
nances are the foundation of Eano resonances in output 


of such complex systems. The transverse field effects on 
EIT presented in Eig.2 and Eig.3 are similar to the single 
and double-Eano resonances but tuned by transverse op¬ 
tical field. We conventionally observe Eano line shapes in 
EIT windows by tunning effective detuning of the system. 
The variation in effective detuning of the system brings 
change to the anti-Stokes scatterings which causes the 
shift in EIT window. In following, we demonstrate Eano 
behavior of system output field with respect to different 
parameters. 

Eig.4 shows single-Eano resonances in the absorp¬ 
tion (real) and dispersion (imaginary) profile of output 
probe field in the absence of transverse field coupling 
Tjeff/ti = 0, as a function of normalized probe field 
detuning Apjujrn and normalized effective detuning of 
the system A/cj^. The coupling of intra-cavity optical 
mode with atomic mode is Gaj^m = 0-03 while, the cou¬ 
pling of optical mode with mechanical mode is kept zero 
{Gml^m = 0)- which means, optomechanical system is 
only coupled to the condensate mode trapped inside the 
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FIG. 6: (Color online) The effects of transverse optical field, 
with coupling r\eff/K — 0.03, on single-Fano resonances are 
presented in the absorption (real) and dispersion (imaginary) 
profiles of output probe field as a function of normalized probe 
field detuning Ap/cjrn and under the influence of normalized 
effective detuning of the optomechanical system A/cjrn- The 
coupling strengths are same as in Fig. 4. (a) shows absorp¬ 
tion quadrature of output optical mode with different detun¬ 
ing strengths Ajum =0.84 (blue curve), A/ujm =0.87 (red 
curve) and Ajum = 0.9 (green curve). While, (b) accounts 
for imaginary quadrature of output field having effective de¬ 
tuning strengths A/uOm — 0.99 (blue curve), A/ujm — 1-02 
(red curve) and Ajum — 1.06 (green curve). Remaining pa¬ 
rameters are same as in Fig. 2. 


cavity. Fig. 4(a) and Fig. 4(b) describe absorption and 
dispersion profile,respectively, in output probe field as a 
function of normalized probe detuning. Blue curve in 
absorption shows Fano line with effective system detun¬ 
ing Ajujm = 0.84. While, red and green curves in real 
quadrature represent fano behavior under influence of ef¬ 
fective detuning A/ujm = 0.87 and A/ujm = 0.9, respec¬ 
tively. Similarly, blue curve in dispersion profile shows 
the existence of Fano resonance with effective system de¬ 
tuning A/ujm = 0.99. Besides, red and green curves in 
imaginary quadrature of output field represent fano be¬ 
havior under influence of effective detuning A/ujm = 1.02 
and A/uJrn = 1-06, respectively. We can observe, each 
curve with different height follow a same dip in absorp¬ 
tion and dispersion response which causes the formation 
of resonance in out-going optical mode. By analyzing 
these curves, one can predict the formation of Fano res¬ 
onance in the output field. 

We further investigate the existence of double-Fano 
resonances in output probe field by introducing another 



FIG. 7: (Golor online) The behavior of double-Fano reso¬ 
nances are demonstrated in the absorption (real) and dis¬ 
persion (imaginary) quadratures of output probe field as a 
function of normalized probe field detuning Aplum and effec¬ 
tive detuning Ajum under the notable effects of transverse 
optical field having strength T]eff/i^ = 0.03. The system 
coupling strengths are same as in Fig. 5. (a) shows absorp¬ 
tion profile of output field under effective detuning strengths 
Ajujm = 0.83 (blue curve), Ajum = 0.85 (red curve) and 
A/ujm — 0.97 (green curve). On the other hand, (b) presents 
dispersion quadrature of output field with effective detuning 
values Ajujm — 1.05 (blue curve), Ajum = 1.1 (red curve) 
and A/ujm — 1.15 (green curve). All other parameters are 
same as in Fig. 2. 


coupling in the optomechanical system and modifing ef¬ 
fective detuning [45]. As the phenomenon of EIT is 
very sensitive to the coupling with different degrees of 
freedom in the system. Therefore by introducing an¬ 
other coupling, we can convert single-Fano resonance to 
double-Fano resonance. Fig. 5 shows such double-Fano 
resonances in the absorption and dispersion profile of out¬ 
put probe field in the absence of transverse field coupling 
Tjeffln = 0 and as a function of normalized probe field 
detuning Apjoj^ and normalized effective detuning of the 
system A/oOm- The coupling of intra-cavity optical mode 
with moving-end mirror is Gat^m = 0.1 and the coupling 
of optical mode with condensate mode is Gat^m — 0.08. 
Fig. 5 (a) describes absorption and Fig. 5(b) describes dis¬ 
persion profile in output probe field as a function of nor¬ 
malized probe detuning. Blue curves, in Fig. 5(a) and 
Fig. 5(b), show the double-Fano line with effective detun¬ 
ing values Ajura = 0.85 and Ajuj^ = 1.1, respectively. 
Similarly, red curves, in absorption and dispersion, ac¬ 
commodate the double-Fano response under the influence 
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of effective detuning /^jb^rn =0.94 and A/cj^ = 0.96, re¬ 
spectively and green curves accounts for the influence of 
effective detuning A/cj^ = 0.87 and A/cj^ = 1.15, re¬ 
spectively. By analyzing these results, we come to know 
the existence of single-Fano resonances as well as double- 
Fano resonances in the output probe field of the system 
[44, 45]. 

In previous Fano resonance results, we have ignored 
the effects of transverse optical field coupling. But it 
will be very important to keep these effects and ana¬ 
lyze the behavior of Fano resonances. Fig. 6 illustrates 
such effects on single-Fano resonances emerging in out¬ 
put field spectra in the presence of transverse field cou¬ 
pling r]eff/i^ = 0.03. Fig. 6(a) shows real quadrature of 
out-going mode where, blue, red and green curves corre¬ 
sponds to the effective detuning strengths A/oOrn = 0.85, 
A/uJm = 0.87 and A/oUm = 0.9, respectively. On the 
other hand. Fig. 6(b) represents imaginary quadrature 
of output field where, blue, red and green curves cor¬ 
respond to the influence of system detuning strengths 
Ajuj^ = 0.99, A/uJm = 1-02 and Ajuj^ = 1.06, respec¬ 
tively. One can observe, how quadratures of singe-Fano 
lines are increased due to the presence of transverse field. 
Transverse optical field causes scattering of photon inside 
the cavity which gives rise to intra-cavity photon num¬ 
ber and this nonlinear factor brings modification to the 
out-going optical mode of cavity. It is understood that if 
we further increase the strength of transverse coupling, 
it will definitely modify Fano behavior in output field. 

Fig. 7 shows similar behavior for double-Fano reso¬ 
nances in real and imaginary profile of out-going probe 
field under the influence of transverse field strength. 
Fig. 7(a) shows absorption and Fig. 7(b) shows dispersion 
behavior at transverse optical field strength = 

0.03. Blue curves, in Fig. 7(a) and Fig. 7(b), show the ef¬ 
fects of r]ef f /n on double-Fano curves appearing in out¬ 
going mode with effective detuning values A/ujm = 0.85 
and A/uJm = 1.1, respectively. While, red curves, in 
absorption and dispersion quadratures, illustrates the 
effects on double-Fano response with effective detun¬ 
ing A/ujm = 0.94 and Ajuj^ = 0.96, respectively and 
green curves accounts for the similar response double- 
Fano resonance under the influence of effective detuning 
Ajuj^ = 0.87 and A/ujm = 1.15, respectively. By com¬ 
paring results of Fig. 4 and Fig. 5 with Fig. 6 and Fig. 7, we 
can easily note the effects of transverse optical field on 
the double-Fano resonance of the optomechanical system. 
The absorption and dispersion quadratures of single-Fano 
as well as double-Fano resonances are notably modified 
by increasing transverse optical field strength which we 
can observe in presented results. 


VI. CONCLUSION 

In conclusion, we discuss the controllability of elec- 
tromagnetically induced transparency (FIT) and Fano 
Resonances in hybrid BEC-optomechanical system which 
is composed of cigar-shaped Bose-Einstein condensate 
(BEC) trapped inside high-finesse Eabry-Perot cavity 
driven by a single mode optical field along the cavity 
axis and a transverse pump field. As, the transverse op¬ 
tical field directly interacts with condensate mode which 
causes the scattering of transverse photon inside the cav¬ 
ity, so, by varying transverse field, we can modify the dy¬ 
namics of system. We have shown the controlled behav¬ 
ior of electromagnetically induced transparency (EIT) in 
output probe field by using transverse field. We dis¬ 
cuss existence of single-EIT window in output field of 
cavity-optomechanics in the absence of moving-end mir¬ 
ror, which means intra-cavity optical mode was only cou¬ 
pled to atomic mode (BEC) of the system. The single- 
EIT as well as double-EIT windows in output probe laser 
field are efficiently amplified by increasing the strength 
of transverse optical field. Furthermore, the single and 
double-Fano resonances are discussed in out-going probe 
field of the system. The influence of transverse optical 
field is also been studied on Fano resonances of the sys¬ 
tem. Moreover, we have suggested a certain set of ex¬ 
perimental parameters to observe these phenomena in 
laboratory. 

In future, we will apply this method of control to dis¬ 
cuss the nonlinear dynamics of hybrid system, especially 
to explore dynamics of interacting Bose-Einstein conden¬ 
sate in such complex systems. Besides, we intend to ex¬ 
tend this method to study the controllability of novel 
phenomenon like entanglement. Additional goals include 
the study of effects of spin-orbit coupling [7, 8] using 
magnetic field in hybrid BEC-optomechanical systems. 


Acknowledgments 

This work was supported by the NKBRSFC under 
grants Nos. 2011CB921502, 2012CB821305, NSEC un¬ 
der grants Nos. 61227902, 61378017, 11434015, SKLQO- 
QOD under grants No. KF201403, SPRPCAS under 
grants No. XDB01020300. We strongly acknowledge fi¬ 
nancial support from CAS-TWAS President’s PhD fel¬ 
lowship programme (2014). 


[1] T. J. Kippenberg and K. J. Vahala, Science 321, 1172 
(2008). 

[2] F. Brennecke, S. Ritter, T. Donner, and T. Esslinger, 
Science 322, 235 (2008). 


[3] H. Ritsch, P. Domokos, F. Brennecke, and T. Esslinger, 
Rev. Mod. Phys. 85, 553 (2013). 

[4] Ying Dong, Jinwn Ye, and Han Pn, Phys. Rev. A 83, 
031608(R) (2011). 



10 


[5] Dae-II Choi and Qian Niu, Phys. Rev. Lett. 82, 2022- 
2025 (1999); W. M. Liu, B. Wu, and Qian Niu, Phys 
Rev. Lett. 84, 2294-2297 (2000). 

[6] W. Vincent Liu, Phys. Rev. Lett. 79, 4056 (1997); M. 
Lewenstein and W. V. Liu, a News V Views article, Nat. 
Phys. 7, 101 (2011). 

[7] Zi Cai, Xiangfa Zhou, and Congjun Wu, Phys. Rev. A 
85, 061605(R) (2012); Zi Cai, Yu Wang, and Congjun 
Wu, Phys. Rev. B 86, 060517(R) (2012). 

[8] L. Dong, L. Zhou, B. Wu, R. Balasubramanian, and Han 
Pu, Phys. Rev. A89, 011602(R)(2014); Hui Hu, B. Ra- 
machandhran, Han Pu, and Xia-Ji Liu, Phys. Rev. Lett. 
108, 010402 (2012). 

[9] V. B. Braginsky, Measurement of Weak Forces in Physics 
Experiments, University of Chicago Press, Chicago 
(1977). 

[10] S. Mancini, V. Giovannetti, D. Vitali, and P. Tombesi, 
Phys. Rev. Lett. 88, 120401 (2002). 

[11] L. Tian and P. Zoller, Phys. Rev. Lett. 93, 266403 (2004). 

[12] A. Naik, et a/.. Nature (London) 443, 193 (2006). 

[13] A. D. O'Connell et al., Nature 464, 697 (2010). 

[14] J. D. Teufel et a/.. Nature 475, 359 (2011). 

[15] J. Chan et a/.. Nature 478, 89 (2011). 

[16] K. Liu, L. Tan, C. H. Lv, and W. M. Liu, Phys. Rev A 
83, 063840 (2011). 

[17] Q. Sun, X. H. Hu, A. C. Ji, and W. M. Liu, Phys. Rev A 
83, 043606 (2011); Yu Yi-Xiang, Jinwu Ye, and W. M. 
Liu, Nature, Scientific Reports 3, 3476 (2013). 

[18] S. Groeblacher, K. Hammerer, M.R. Vanner, and M. As- 
pelmeyer. Nature (London) 460, 724 (2009). 

[19] J. D. Teufel, D. Li, M. S. Allman, K. Cicak, A. J. Sirois, 
J. D. Whittaker, and R. W. Simmonds, Nature (London) 
471, 204 (2011). 

[20] E. Verhagen, S. Deleglise, S. Weis, A. Schliesser, and T. 
J. Kippenberg, Nature (London) 482, 63 (2012). Dae- 
II Choi and Qian Niu, Phys. Rev. Lett. 82, 2022-2025 
(1999). 

[21] K. Zhang, W. Chen, M. Bhattacharya, and P. Meystre, 
Phys. Rev. A 81, 013802 (2010). 

[22] Ying-Dan Wang and A. A. Clerk, Phys. Rev. Lett. 108, 
153603, (2012). 

[23] S. Singh, H. Jing, E. M. Wright, and P. Meystre, Phys. 
Rev. A 86, 021801, (2012). 

[24] Gabriele De Chiara, Mauro Paternostro, and G. Massimo 
Palma, Phys. Rev. A 83, 052324 (2011). 

[25] M. Abdi, S. Pirandola, P. Tombesi, and D. Vitali, Phys. 
Rev. Lett. 109, 143601 (2012). 

[26] M. Abdi, S. Pirandola, P. Tombesi, and D. Vitali, Phys. 
Rev. A 89, 022331 (2014). 

[27] Eyob A. Sete, H. Eleuch, and C. H. Raymond Ooi, J. 
Opt. Soc. Am. B 31, 2821 (2014). 

[28] T. Shi, Longhua Jiang, and Jinwu Ye, Phys. Rev. B 81, 
235402 (2010); A. M. Dudarev, R. B. Diener, B. Wu, 
M. G. Raizen, and Q. Niu, Phys. Rev. Lett 91, 010402 


(2003). 

[29] K. A. Yasir, M. Ayub, and E. Saif, J. Mod. Opt. 61, 1318 
(2014); M. Ayub, K. A. Yasir, and E. Saif, Laser Phys. 
24, 115503 (2014). 

[30] A. C. Ji, X. C. Xie, and W. M. Liu, Phys. Rev. Lett. 99, 
183602 (2007); A. C. Ji, Q. Sun, X. C. Xie, and W. M. 
Liu, Phys. Rev. Lett. 102 , 023602 (2009); L. Tan, Y. Q. 
Zhang, and W. M. Liu, Phys. Rev. A 84, 063816 (2011). 

[31] Marian O. Scully and M. Suhail Zubairy, Quantum Op¬ 
tic, Cambridge University Press (1997). 

[32] A. H.Safavi-Naeini, et a/.. Nature 472, 69 (2011). 

[33] S. E. Harris, Phys. Today 50, 36 (1997). 

[34] S. E. Harris, J. E. Eield, and A. Imamoglu, Phys. Rev. 
Lett. 64, 1107 (1990). K. J. Boiler, A. Imamoglu, and S. 
E. Harris, Phys. Rev. Lett. 66, 2593 (1991). 

[35] M. M. Kash, et al, Phys. Rev. Lett. 82, 5229 (1999); M. 
D. Lukin and A. Imamoglu, Nature 413, 273 (2001). 

[36] G. S. Agarwal and S. Huang, Phys. Rev. A 81, 041803(R) 

( 2010 ). 

[37] M. Asjad, J. Russ. Laser. Res. 34, 3 (2013); M. Asjad, J. 
Russ. Laser. Res. 34, 2 (2013). 

[38] D. E. Chang, A. H. Safavi-Naeini, M. Hafezi, and O. 
Painter, New J. Phys. 13, 013017 (2010); A. H. Safavi- 
Naeini and O. Painter, New J. Phys. 13, 013017 (2011). 

[39] S. Weis, et al, Science 330, 1520 (2010). 

[40] E. Massel, et a/., Nat. Commun. 3, 987 (2012). 

[41] U. Eano, Phys. Rev. 124, 1866 (1961). 

[42] B. H. Bransden and C. C. Jean Joachain, Physics of 
Atoms and Molecules, 2nd ed. (Addison-Wesley, New 
York, 2003), Chap. 4. 

[43] Benjamin Gallinet and Olivier J. E. Martin, Phys. Rev. 
B 83, 235427 (2011). 

[44] K. Qu and G. S. Agarwal, Phys. Rev. A 87, 063813 
(2013). 

[45] M. J. Akram, E. Ghafoor, and E. Sair, J. Phys. B, 48, 
065502 (2015). 

[46] J. Estve, C. Gross, A. Weller, S. Giovanazzi, and M. K. 
Oberthaler, Nature (London) 455, 1216 (2008). 

[47] E. Brennecke, et al., Nature (London) 450, 268 (2007). 

[48] Mauro Paternostro, Gabriele De Chiara, and G. Massimo 
Palma, Phys. Rev. Lett. 104, 243602 (2010). 

[49] J. M. Zhang, E. C. Cui, D. L. Zhou, and W. M. Liu, 
Phys. Rev. A 79, 033401 (2009). 

[50] Kashif Ammar Yasir, and Wu-Ming Liu, 
arXiv: 1503.08752 [quant-ph] (2015). 

[51] Shuai Yang, M. Al-Amri, Jrg Evers, and M. Suhail 
Zubairy Phys. Rev. A 83, 053821 (2011). 

[52] X. E. Zhang, et aL, Phys. Rev. Lett. 110, 090402 (2013). 

[53] C. K. Law, Phys. Rev. A 51, 2537 (1995). 

[54] C. W. Gardiner, Quantum Noise (Berlin, Springer, 1991). 

[55] M. Paternostro, et aL, New Journal of Physics 8, 107 
(2006); V. Giovannetti and D. Vitali, Phys. Rev. A 63, 
023812 (2001). 


